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Abstract. This paper studies a deflation method to reduce and to 
solve linear differential-algebraic equations (DAEs) . It consists to define 
a sequence of DAEs with index reduction of one unit by step. This is si- 
multaneously performed by substitution and differentiation. At the end 
of process, we obtain at most an ODE and a list of algebraic constraints 
which solve the initial DAE. We show on classical examples how works 
the method. Moreover, we explain how this method extends in the case 
of linear time- varying DAEs. 
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1. Introduction 

This paper focuses on the study of a deflation type method to solve lin- 
ear differential-algebraic equations (DAEs). Deflation type methods arise 
in many problems involving eigenvalues, roots of polynomial systems etc... 
Roughly speaking, a deflation process replaces the original problem by a 
problem of smaller size. In the case of DAEs, it consists to determine a 
sequence of DAEs of size strictly decreasing ; at the end of this process, we 
obtain at most an ODE and a list of algebraic constraints satisfied by the 
general solution of the DAE. In this paper, we deal with both linear time- 
invariant DAEs and linear time- varying DAEs. 

The literature on linear DAEs is extremely rich with a large classical docu- 
mentation. In the recent past, the monographs of Kunkel and Mehrmann [9] 
and Riaza [17] give an excellent overview and permit to understand the dif- 
ferent approaches from the works of pioneers, Weierstrass and Kronecker. 
The linear case is related to properties of matrix pencils, which are well de- 
scribed in [2], where the notion of regularity provides fundamental results. 
The diversity of perspectives and techniques for the DAEs is probably due 
to the central notion of index. In fact, if the notion of Kronecker index in the 
linear time-invariant DAEs' case is an unifying concept, it is not the same 
for linear time-dependent case. Several index notions appear according to 
the point of view adopted : differentiation index developed by Campbell, 
Gear and Petzold [1] ; reduction methods and geometric index studied by 
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Rheinboldt and Rabier [16] and [8] ; projector-based methods and tractabil- 
ity index introduced by Griepentrog and Marz [5] ; strangeness index by 
Kunkel and Mehrmann [9] ; perturbation index by Hairer [8] ; structural in- 
dex by Jacobi [12] and Pryce [13]. Except for the structural and strangeness 
indices, all these indices are equal in the linear time-invariant DAEs' case. 
The complexity of the deflation method, which is proposed here, deals with 
the Kronecker index (time-invariant case) and the differentiation index (time- 
varying case) . The technical background of the method uses elementary al- 
gebra as LU decomposition or singular value decomposition performed on 
the matrix E of the following DAE : Ex = Ax + f. 

In order to provide some motivations at our work, we briefly remember the 
main contributions in the field of linear DAEs with constant coefficients. 
The first remark concerns the regularity assumption for linear DAE. In fact, 
if we don't use this assumption, there are more than countably many dif- 
ferent solutions as it is well described in [5] pages 14-15. The classical re- 
duction of this DAEs type deals with the Weierstrass canonical form and 
Kronecker index . Then the linear DAE is decoupled in two subproblems 
of type u = Ju + g and Nv = v + h, where the matrix J and the nilpo- 
tent matrix N are related to the Jordan form of (XE + A)~ 1 E, for some 
A € K. From a computer algebra point of view, this way permits to obtain 
an expression of the solution, i.e v 

(*) = - T^o Nk h (k) {t) where v is the 
Kronecker index, i.e N u ~ l ^ and N v = 0. But from a numerical analysis 
point of view, finding the Jordan normal form is known to be numerically 
unstable : small errors can make large differences in the result, see [4] or 
the discussion in [3] sec. 7.6.5. We also note the Weierstrass canonical form 
needs to find A such that XE + A is invertible and to compute (XE + A)~ 1 E 
before to perform a Jordan normal form. These numerical drawbacks led us 
to propose a reduction of type deflation. 

On the other hand, the differentiation method, described in [1] page 20, is 
based on the differentiation of constraints. A new DAE is obtained and the 
procedure is repeated until to obtain an ODE. The result (Theorem 2.3.3 [1]) 
says that the number of steps is precisely the Kronecker index in the time- 
invariant case, and each differentiation of constraints reduces the index by 
one. But the final ODE has many additional solutions. Instead the deflation 
method proceeds by substitution, and the differentiation is performed during 
this stage. We will state that the number of steps is the minimum of the 
Kronecker index and the rank. Moreover, solutions of the reduced system 
are those of the initial system. There are other types of reduction in the 
literature ; mainly the one introduced in [5], based on projector methods 
and generalized in [6] for DAEs with higher index. A good introduction of 
these projector methods is done in [17] : the solution of the DAE is expressed 
thanks to matrix chains, which length is equal to the Kronecker index. How- 
ever, our method appears to be technically simpler. 

More recently in [18], an index reduction based on substitution method has 
been developed for DAEs of the type Ex = Ax + f(t) with 
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where B is an invertible matrix. Introducing the Schur complement D = 
M — LB~ X K, we obtain a new DAE 

Dx 2 = h~ LB- l h (1) 

under the constraints 

Bx x = h- K± 2 . (2) 
It is proved the index of the DAE 1 is one lower than that of initial DAE. 
Next the DAE 1 is numerically solved and the solution of the initial DAE is 
recovered thanks to the system of linear equations 2. This method is mainly 
proposed in the context of electrical circuit where it is possible to find a non- 
singular constant submatrix B of XE + A. Certainly a reduction of an index 
unit improves the numerical results, but the drawback of this method is the 
DAE with higher index. In some sense, our deflation method generalizes this 
approach. Finally, we have to mention the paper of Linh and Mehrmann [10] 
where a similar transformation of the DAE and the Schur complement are 
used in the strangeness- free context. 

In this paper, we will separately treat the time-invariant case and the time- 
varying case because of their structural differences. In section 2, we present 
the deflation method in the time-invariant case and the main results, which 
will be proved in section 3. The section 4 gives some examples which il- 
lustrate the method. Finally the section 5 shows how works the deflation 
method in the time-varying case. 

2. Deflation method for linear time-invariant DAEs. 

A linear time-invariant DAE is a problem of the form : 

Ex(t) = Ax(t) + f(t), (3) 

where E and A are constant matrices in M nxn and f : I —■ W 1 is a sufficiently 
smooth function, defined on an open interval ICR. We assume the rank 
of E is r and this DAE is regular : there exists A € R such that the matrix 
XE + A is non-singular [2] . 

Using both basic linear algebra and substitution of certain variables, the 
main idea is to separate the differential part and the algebraic part of the 
DAE 3. This goal can be attained using a decomposition of the matrix as LU 
decomposition or singular value decomposition (SVD) which are well studied 
in [3]. For example if the LU decomposition is used, E = LU and the DAE 3 
is equivalent to 

Ux(t) = L~ x Ax(t) + L- X f{t). 
Hence the algebraic constraints appear since there are n — r zero rows in 
the matrix U. A similar result is obtained using SVD. For this reason, we 
present the deflation method without reference to the decomposition used 
to factorize the matrix E. The following lemma is the key of the deflation 
method. 

Lemma 1. There exist an invertible matrix U and a permutation matrix P 
such that XE + A reads 

XE + A = U ( XS + K AT + L V> 
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where S,K G M rxr ; T,L G Rrxn-r. M G M"-^ and JV G 
non-singular matrix. 

Proof. Since the rank of E is r , we can factorize E = UT, := U I ^ ) where 



f7 is invertible and .F of size r x n with rank r. Writing U 1 A 
where Ai G W xn and A 2 G R™- rx ™, we obtain 

XF + Ai 





^2 



A£ + A = U , , 

^2 

The matrix ^2 has full rank since \E + A is invertible. Hence there exists 
an invertible submatrix N G ]J ra - r ' xri - r G f A2. Introducing an appropriate 
permutation matrix P, we obtain the desired formula. n 

With the notations above, it is easy to see the DAE 3 is equivalent to 

where U~ 1 f = ^ ^ V After an easy computation, the DAE 3 is finally 
reduced to 

(S - rJV -1 M)ti = (-K" - LiV -1 M)u + TA^A - LN~ l h + 5 (5) 
v = -N-i-Mu - N^h (6) 
P- 1 ^ = (-u,t;) r (7) 

The size of the DAE 5 is less than that of the DAE 3. This conduces to 
define a deflated DAE of an initial DAE. 

Definition 1. We note 

Ei = S- TN~ l M, 



A X = K - LN~ X M, 

f 1 = TN^h - LN^h + g, 

x 1 = u. 

We say that Eix 1 = A\x l + f\ is a deflated DAE of 3. Moreover P is 
the permutation matrix and Mu + Nv + h = is the algebraic constraint 
associated to this deflation. 

Let us remark there are in general several deflated DAEs of a given DAE : 
in fact it depends upon the choice of matrices M and N . Fortunately the 
ranks of matrices E\ and A\ are invariant as it is stated in the following. 

Theorem 1. Let us consider two deflated DAEs of the DAE 3 represented 
respectively by the matrix pencils (E\,Ai) and (E\,Ai). Then we have 
rank(Ei) = rank(Ei) and rank(A\) = rank(Ai). 

Moreover, the regularity is preserved during a step of deflation. 

Theorem 2. If \E + A is non-singular then \E\ + A\ is non- singular. 
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From the two previous theorems which will be proved in the next section, 
it follows first that the deflation process is independent of the choice of the 
two matrices M and N, and next that the regularity is preserved. This sug- 
gests the following reduction procedure. 

Deflation algorithm. 

Input : E = E, A = A, f = f, x° = x, r_i = n. 

Step j + l,j> 0. 

If Ej is singular and Ej ^ 

1- Let rj the rank of Ej . 

2- Compute (Ej+x, Aj+i, fj+i) from (Ej,Aj,fj) using the formulas of 
the definition 1. Let Pj the permutation matrix of this reduction. 

3- Compute the change of variable (x\,X2) T ■= Pf x ^ where the size 
of x\ is rj. 

4- Compute the algebraic constraint : = Mjx\ + NjX J 2 + fj,rj+i-.r j ^ 1 , 
where /j r +i r _i means the coordinates r,- + 1 to Tj-x of the vector 
fj. 

5- Let x 3+ := x\. 
else stop. 

Output For j > the sequences 

DAEs algebraic constraints 

E j+1 xi +1 = A j+1 xi +1 + f j+1 , = M jX { + N jX { + /,.,,;. , . 

The main goal of this paper is to prove the deflation algorithm stops in a 
finite number of steps. More precisely 

Theorem 3. The number of steps of the deflation algorithm is bounded by 
min(rank(E) , ind(E , A)) where ind(E, A) is the Kronecker index of the pen- 
cil (E, A) . Moreover, the ranks of matrices Ej and Aj determine a sequence 
of invariants which are characteristic for E and A. 

The general solution is described by the following result. 

Theorem 4. Let k be the number of steps of the deflation algorithm. Then 
the coordinates of the solution satisfy 

E k x k = A k x k + f k (8) 

where E^ is invertible or equal to zero and 

x{ = -N^Mjxi - Nf 1 f j , rj+1:rj _ 1 , < j < k - 1. (9) 

Proof. From the definition of the deflation algorithm. 

3. Proofs of theorems of the section 1. 

To prove that let us remember some fundamental notions. The index of 
a matrix B is the smallest integer which verifies Ker(B k ) = Ker(B k+1 ). If 
the index B is equal to k then it is equivalent to the transversality condition 
or the range-nullspace decomposition [11] page 394: 

R n = Ker{B k ) ®Im{B k ). 
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If the index of B is zero then B is non-singular. The Kronecker index of the 
DAE 3 is the index of the matrix {\E + A^E. We denote by ind(E,A) 
this index. Moreover, lemma 7 page 196 of [5] shows that ind(E, A) = 
ind(EQ, AQ), if Q is invertible. 
We first prove the theorem 1. 

Proof of theorem 1. Let us suppose there exist two choices of matri- 
ces (S,T,K,L,M,N) and (S,f,K,L,M,N) for which the lemma 1 holds. 
Then there exists a permutation matrix Pq such that 

S T\ ( S f 
M N J \ M N 

Since N and N are non-singular, we can write 

TN- 1 I \ ( N \ / N- l M I \ 
I ) { Eh ) V / J Po 

TN- 1 I \ f N \ / N- l M I 

i o y V o y V 1 

Hence the ranks of the matrices E\ and E\ are equal. By a same way, we 
show rank(Ai) = rank(Ai). □ 
We now prove the theorem 2. 

Proof of theorem 2. From lemma 1 it follows the matrix ASP + C/ _1 AP 
is non-singular. Using the Schur complement of this matrix, we have 



XZP + U->AP=( XS + K XT + L ) 



(XT + LjN- 1 I \ f N \ f N~ X M I 
I ) [ XEi + At ) [ I 



It follows the matrix A-Ei + A\ is non-singular. □ 
To prove the theorem 3 we need some lemmas. 



1- (XEP + U 1 AP) l Y>P - | , | 1/Y ,-|. 



Lemma 2. With the notations of the section 1 and E = UT,, C\ = \E\-\-A\, 
we have : 

-N^MC^S -N- l MC^ L T 

2- 

((AEP +U- 1 APY l HP) k 

{C^E^^C^S {C^E^^C^T 

Proof. A straightforward computation gives the result of the part 1. 
For the part 2, we apply the lemma below to the identity of the part 1. □ 

Lemma 3. Let k > 1. 

(A B \ k _ ( {A + BC) h ~ l A {A + BC) k - l B \ 
\CA CB ) ~ I C(A + BCf^A C(A + BC) k ~ l B ) 
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Proof. By induction. □ 

We next state a result which appears in [5], in theorem 13 page 198. We 
remember the proof for sake of completion. 

Lemma 4. The index of the DAE 3 is equal to one iff x £ Ker E and 
Ax G Im E =>■ x = 0. 

Proof. Let us suppose the index of the DAE 3 is one. This implies 
M. n = Im(XE + A)~ l E Ker(XE + A)^E. Let us suppose x G Ker E and 
Ax G ImE. Then x G Ker E implies x G Ker (XE + A)~~ 1 E. On the other 
hand, there exists y G 1" such that Ax = Ey. Since (XE + A)x = Ax and 
XE + A is non-singular we have x = (XE + A)~ 1 Ey G Im (XE + A)~ 1 E. 
Hence x G Ker(XE + A)" 1 ^ n Im(XE + A)" 1 ^. It follows x = 0. 
Conversely, let us suppose x G Ker E and ^4x G ImE => x = 0. Assuming 
the index A; > 1, there exists y ^ such that (A-E 1 + ^4) _1 £'y ^ and 
((AE + A)" 1 ^) 2 y = 0. Hence z := (XE + A)~ 1 Ey G Ker.E. Prove that 
Az G ImE. In fact (Ai? + = hence Az = Ey G ImE. This means 
z = in contradiction with the definition of z. □ 

Lemma 5. The index of the DAE 3 is equal to one iff E\ is non- singular : 

Proof. Let us suppose ind(E,A) = 1 and prove that KerE\ = {0}. 

( S T \ 

It is equivalent to show that the matrix I I is non-singular since 



v M N / 

the Schur complement of this matrix is E\. Let us consider x such that 
S T \ 

M N ) X = ^' ^ i m P nes both x G Ker(S, T) = Ker(T,P) and (M, N)x = 
0. Hence ( £ M x _ ( <* £ >* ). Since the rank of EP is ,, there 
exists y G 1" such that f (K,L)x \ _ yjp^ j n £ ac ^. we ^ave x g Ker(T,P) 
and U~ 1 APx = T,Py G Im(T,P). From lemma 4 it follows x = and the 
matrix f ^ ^ ~\ is non-singular. 
Let us suppose now -Ei is non-singular. 

From lemma 1, we have ind(E,A) = ind(YiP,U~ l AP). From lemma 4, 
proving that ind(EP, U~ 1 AP) = 1 is equivalent to establish the assertion 
x G Ker(Y,P) and U^APx G Jm(SP) =>- x = 0. If we have x G Ker(T,P) 
and C/^APx G Jm(SP), this implies (5 , ,T)x = and (M,N)x = 0. Since 
5 T 

is non-singular, it follows x = 0. □ 



M N 

We can now state the result which links the indices of the DAE 3 and the 
deflated DAE 5. 

Theorem 5. If E x / then ind(E u A x ) = ind(E,A) - 1. 

Proof. The case ind(E, A) = 1 is treated by the lemma 5. 

Let us suppose first that k = ind(E, A) > 2. Since (XE+A)~ l E = P (XT.P+ 

U- 1 AP)- 1 Y.PP- 1 , the indices ind(E,A) and ind(T,P, U^AP) are equal. 
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We then have 

R n = Ker ((ASP + U^APy^pf Im ((ASP + U^APy^pf . 

(10) 

It is sufficient to show 

x G Ker {C^Exf- 1 n Im (Cf 1 ^)* -1 => a; = 0, 

where d = A#i + A x . If x € Ker (Cf 1 ^)^ 1 then (Cf 1 ^)*- 2 ^- 1 ^ - 
TN~ 1 Mx) = 0. From the identity of lemma 2 part 2, it follows 

(x, -N~ l Mx) T G Ker ((ASP + U~ 1 APy 1 Y l P) k . 

Now if x € Im (Cf 1 Pi) fc " 1 i there exits y G M r such that x = (C^E^^y. 
Since the rank of SP is equal to r, there exists (u, v) T G M. r x R n_r such 
that y = C^iSu + Tv). Then x = (C 1 _1 Pi) fc - 1 Cf 1 (Su + Tu). Always from 
lemma 2 part 2, it follows 

(x, -N" l Mx) T G Im ((ASP + U^APy^pf . 

From 10 we deduce (x, —N~ l Mx) T = and finally x = 0. □ 
Next, we state some properties concerning the rank of E\. 

Lemma 6. Let us consider the regular DAE 3. 

1- Ifind(E,A) > 1 then rank(E) > rank{E\). 

2- Ifind(E,A) = 1 then rank(E) = rank{E\). 

3- Ifind(E,A) > 1 and rank(E) = 1 then E 1 = 0. 

Proof. Prove the part 1. From construction r = rank(E) > rank(Ei) 
since E\ G M rxr . If rank{E\) = r, this implies E\ is non-singular and 
we have ind(E, A) = 1 from lemma 5. This contradicts the hypothesis 
ind(E, A) > 1. The part 2 is a direct consequence of part 1. For the part 3, 
we observe that E\ G M. since rank(E) = 1. If E\ ^ then E^ 1 exists. This 
contradicts ind(E, A) > 1. Hence E\ = 0. □ 
Now, we are able to prove the theorem 3. 

Proof of theorem 3. At each step of the algorithm, the index is strictly 
decreasing from theorem 5. In the same way, from lemma 6, the rank of each 
Ej is also strictly decreasing. The reduction algorithm stops when the index 
or the rank is 0. The result follows. □ 

4. Examples 

We present in the examples below, the different steps of the deflation al- 
gorithm with the differential part in the left and, the algebraic constraints 
in the right. 

Example 1. [1] page 19. 

1 \ / -1 \ 
1 \x= -1 x + f(t). 
000/ \ -1 / 

The Kronecker index is 2. The sequence of DAEs and the algebraic con- 
straints given by the deflation algorithm is successively described in the two 
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steps below. 



Step 1. The permutation matrix consist to swap the columns 2 and 3. Hence 
, I Q\ ( 



S - TN~ M - , 
and the constraints are : 



1 



-1) ( ) = S. The deflated DAE 



1 




xi 
x 3 



-1 



x + 



•ft . 

h - h 



= X 2 -f 3 



Step 2. 



-S3 + h~ h + h 



= xi - / 2 + / 3 



In this case the algorithm directly gives algebraic constraints which deter- 
mine the solution. 

Example 2. N^x = x where N is an elementary nilpotent matrix of size 
/0 1 \ 

kxk : N k = ' ' ■ ' ' • .We have iV^ 1 / and N k = 0. 

1 

V 0/ 

The deflated DAE is Nk_iX\±_i = X\±_i and the constraint relation x k = 0. 

Example 3. This class of DAEs appears in [15] and has a Kronecker index 
equal to 1. 



X2 + x 3 
X2 + X 3 

x 4 + x 5 
x 4 + x 5 




-xx + flit) 

~X2 + f2(t) 
-X3 + f 3 {t) 
-x 4 + fi(t) 

-x 5 + h{t) 



The deflation algorithm stops after one step. 
Step 1. 



X 2 +X A = -X2 + h - h + U q 
±4 = ~Xa + /4-/5 




h 

h + h 

/1 + /2 



Example 4. This is the example 7 [18] which described an electric circuit 
with index three provided by [7]. This DAE is first reduced to a DAE with 
index two and next the solution is numerically computed thanks to the DAE 
solver in Matlab. For this example, our method finds an ODE which gives 
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the exact solution. In our context this DAE reads 

= —x\ — X2 
= —X5 + xq 
= -x 4 + V(t) 

Cx 8 = X 2 

Lx 7 = xq 
= — ax\ — x% 
x 3 + x 7 = 
X4 — ig = 

The index of this DAE is equal to 3. The sequence of deflated DAEs given 
by the algorithm is given in the three steps below. 



Step 1. 



Step 2. 



xj — axi = 
-i 8 = ~V(t) 

Lx-j = Xq 
CX8 = —X\ 
Lx-j = Xq 

-x 8 = -V(t) 
x 7 = -aCV{t) 
= xq + LaCV{t) 



Step 3. 



-i 8 = -V(t) 
x 7 = -aCV(t) 

Consequently the solution satisfies 

- i 8 = 



/ X 5 - Xq \ 

Xi + X 2 

ax\ + x 3 

V X 4 - V(t) J 



= Xl + CV(t) 



= xq + LaCV(t) 



-V(t) 



x 7 = -aCV{t) 
= xq + LaCV(t) 
= xi + CV(t) 

= X 5 - Xq 
= Xi + X 2 

= ax\ + X3 
= x A - V{t) 



5. Deflation method for linear time- varying DAEs 

Now, we consider a linear time-varying DAE 

E(t)x(t) = A(t)x(t) + f{t) (11) 

where E(t), A(t), f(t) are matrices (resp. vector) in M nxn (resp. IR n ) , 
sufficiently smooth, defined on an open real interval /. We assume the rank 
of E(t) is equal to r on the interval /. The assumption of regularity of the 
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DAE 11 on the interval / denned in the section 1 does not apply in this case. 
The classical example [9] page 56 

shows there exists a solution of this equation whereas the determinant of 
XE + A is zero for all A. For this reason, we precise the notion of regularity 
that we will use. It has been introduced under another form in definition 3.1 
of [14]. 

Definition 2. The DAE 11 is geometrically regular on the open interval I 
if 

1- the rank of E(t) is constant on I. 

2- the conclusion of lemma 1 holds. More precisely, there exist a permu- 
tation matrix P and matrices S(t),T(t),K(t),L(t),M(t),N(t) suffi- 
ciently smooth on I such that N(t) is invertible and 

The DAE of the example above is geometrically regular on R + . The 
following definition specifies the formulas for one step of the deflation process. 

Definition 3. Let us suppose the DAE 11 is geometrically regular on the 
interval I. With the notations of the section 1 and U(t)~ 1 f(t) = (g(t), h(t)) T 
we define : 

E X = S- TN~ l M 

A\ = K — LN~ X M + TN~ X {M - NN~ X M) 

f x = TN~ l (h - NN^h) - LN^h + g 

P~ l x = (u,v) T , x 1 = u 

where all the matrices and vectors above depend on t, unless the matrix P. 
We say that Ei{t)x l = Ai^x 1 + fi(t) is a deflated DAE of E(t)x = A(t)x + 
f(t). Moreover P is the permutation matrix and Mu + Nv + h = is the 
algebraic constraint of this deflation. 

Unlike the case of linear time-invariant DAEs, the deflated DAE is not 
necessarily geometrically regular. For example if 



h- > 
















h- > 


1 


■ 


, Mt) = | 


i 





i 


V o 











t 





then 

£iW=(° J), A 1 (t) = 0. 

Consequently we need to suppose the regularity of the deflated DAE to 
continue the deflation process. In this way the deflation algorithm defined 
for linear time-invariant DAE is transposable for linear time-varying DAE. 
Then, we have 
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Theorem 6. Let k be the number of steps of the deflation algorithm applied 
to the linear time-varying DAEs' case. Suppose that each DAE Ej(t)x J = 
Aj(t)x J + fj(t), < j < k — 1 obtained in the deflation algorithm is ge- 
ometrically regular on I. Then the differentiation index of the DAE 11 is 
equal to k. Moreover this index is bounded by the rank of E(t) and the ranks 
of matrices Ej(t) and Aj(t) determine a sequence of invariants which are 
characteristic for E(t) and A(t). 

Proof. From the definition of the differentiation index given in [1] page 
22, it follows this index is equal to k. The other properties follow of the 
definition of the deflation process. □ 



Example 5. This example is treated in [17] page 65. We suppose the 
functions Ci(t), C 2 (t) and L(t) do not vanish. 

— (Ci(t)xi) = x 4 - x 5 

^-(C 2 (t)x 2 ) = -x 3 - x 4 
at 

^(L(t)x 3 ) = x 2 

= X\ — X 2 + Rl(t)xi 

= x\ — R 2 (t)x 5 

We note by E and A the matrices of this DAE which the index depends to 
the values of the parameters. 

Index one. Rx(t) ^ 0, R 2 (t) ^ 0. 

Step 1. We permute the columns 1 and 4 and the columns 2 and 5 in the 
matrices E and A simultaneously. 

CiR 2 x 5 = x 4 - (1 + C 1 R 2 + CiR 2 )x h 
C 2 Ri± A + C 2 R 2 x 5 = -x 3 - (1 + C 2 Ri + C 2 Ri)x A - (C 2 R 2 + C 2 R 2 )x 5 
Lx 3 = -Lx 3 + R1X4 + R 2 x 5 



R1X4 + X\ — x 2 
-R 2 x 5 + x\ 



= 

At this step we have x 1 := (x4,Xs,X3) T and x 2 = (x\^x 2i 



Index two. R\{€) = 0, R 2 (t) 7^ 0. Let us also suppose C\ + C 2 7^ 0. 

Step 1. We permute the columns 1 and 4 and the columns 2 and 5 in the 

matrices E and A simultaneously. 

CiR 2 ± 5 = x 4 - (1 + C\R 2 + C x R 2 )x h / _ \ 

C 2 R 2 x 5 = -x 4 - {C 2 R 2 + C 2 R 2 )x 5 - x 3 = f _^ Xi J 

Lx 3 = R 2 x 5 - Lx 3 

We note E\ and A\ the matrices of this DAE and we have x 1 := (x 4 , X5, x 3 ) T 
and x 2 = (xi,X2) T . 
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/ 1 -d/C 2 

Step 2. We note E x := PE X and A x := PA X with P = 1 

\ 1 
We permute the rows 1 and 3 in the new matrices E\ and A\. 

L(l + c)cx4 + Lbcxs = — (L(c + (?) — Lc)x4 — (R2C 2 + Lcb + Lcb — Lcb)xs 
C1R2X5 = X4 + (ca + 6)x5 

= (1 + C)X4 + 6x5 + CX3. 

with c := C1/C2, a := -C 2 P 2 - Q2P2 and 6 := -1 - CiP 2 + C 2 R 2 c. The 
matrix £2 is invertible since we have C\(t) + C 2 (t) / from the assumption. 
Finally the solution satisfies 

L(l + 0)0x4 + Lbcx§ = — (L(c + c 2 ) — Lc)x4 — (R2C 2 + Lcb + Lcb — Lcb)x§ 
CiR 2 x 5 = x 4 + (ca + 6)x 5 

= (1 + C)X4 + 6x5 + CX3 

= x\ — x 2 

= -R2X 5 + x\. 



Index three. R x (t) = 0, R 2 {t) / and d(t) + C 2 (t) = 0. Hence Ci(t) + 
C 2 (t) = 0. 

Step 1. It is the same as step one for index two. 

CiR 2 x 5 = x 4 - (1 + Cii? 2 + CiP 2 )x 5 
C\R 2 x 5 = x 4 - ((7ii?2 + £1^2)^5 - %3 = 
Lx 3 = R 2 x 5 - Lx 3 

We note E\ and the matrices of this DAE and we have x 1 := (X4, X5, x^) T 
and x 2 := (xi, X2) T . 

Step 2. We proceed as in step 2 for index 2. 

Lx 5 = -(i?2 + L)x 5 
CiR 2 x 5 = x 4 - (1 + Cii? 2 + C 1 R 2 )x 5 
= x 5 + x 3 . 

Now, the matrix E 2 is not invertible. We have ) T and x 2 := X3. 

Siep 5. We note E 2 := PE 2 and A 2 := PM with P = 

We next permute the rows 1 and 2 and columns 1 and 2 in the new matrices 
E 2 and A 2 . Then 



( xi-x 2 \ 



-L/(dP 2 ) 

1 



LCiR 2 x$ = —(La + 6)x5 
= 6x5 + LX4 
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with a = 1 + dR 2 + CiR 2 and b = C X R\ + C X R 2 L -L- dR 2 L - C X R 2 L. 
Finally the solution satisfies 

LC1R2X5 = —(La + b)xFj 
= bx$ + LX4 
= x 5 + x 3 

= X\ — X2 

= -R2X 5 + xi. 



6. Conclusion 

In this paper, a deflation algorithm for DAEs has been studied. It high- 
lights the key role of index notions, especially the Kronecker index and the 
differentiation index. Moreover, the construction of the algorithm establishes 
a connection between the rank of E and the index of the DAE. There is the 
choice to use LU decomposition or SVD at each deflation step. The invari- 
ants of this algorithm are the index of the DAE and the successive ranks T{ 

of deflated DAEs. The arithmetic complexity of the deflated sequence is in 

k 

0( ^2 r f) operations. This reduction provides at most an ODE, which is 

i=-l 

solvable with classical techniques, and also algebraic equations. This algo- 
rithm has the advantage of being technically simple and therefore it brings an 
additional process in order to reduce and to solve linear differential-algebraic 
equations. 
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